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ABSTRACT 


In  acoustics,  the  so-called  Transition,  or  T-matrix  relates  the  incident  and  scattered 
acoustic  pressures  of  an  object  or  scatterer.  The  T-matrix  of  a  thin  steel  spherical  shell  in 
water  has  been  determined  by  the  COMSOL  Multiphysics  Finite-Element  Code.  The 
shell  has  an  outer  radius  of  0.5m  and  a  thickness  of  1cm.  It  is  driven  at  a  frequency  of 
474  Hz  such  that  ka=  1  (where  k  is  the  acoustic  wave  number  and  a  is  the  outer  radius  of 
the  shell).  A  standing  wave,  represented  by  a  spherical  Bessel  function,  is  incident  onto 
the  shell  surface  and  the  corresponding  scattering  coefficient  is  computed.  The  approach 
is  divided  into  three  portions.  Firstly,  a  fluid-loaded  rigid  sphere  is  modeled  using  the 
Acoustic-Shell  Interaction  (ACSH)  physics  mode  to  examine  the  functionality  of 
COMSOL.  It  also  explores  the  degree  of  improvement  when  a  refined  fluid  mesh  is 
applied.  Secondly,  a  thin  spherical  shell  is  modeled  in  the  ACSH  physics  mode.  This  will 
examine  the  credibility  of  COMSOL  to  obtain  accurate  results  based  on  thin  shell 
approximation.  Finally,  a  true  3D  finite  element,  employing  the  3D  elastic  theory,  is 
created  using  the  Acoustic-Structure  Interaction  (ACSI)  physics  mode.  The  resulting 
diagonal  T-matrix  elements  achieved  an  accuracy  of  0.1%  relative  to  the  analytical  T- 
matrix.  Ultimately,  these  results  will  be  applicable  to  the  modeling  of  the  radiation  from 
an  arbitrarily  densely-packed  array  of  sonar  transducers. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  the  purview  of  underwater  naval  operations,  accurate  target  detection  motivates 
the  need  to  deploy  active  sonar  arrays.  Understanding  the  interactions  among  the 
elements  of  an  array  of  sonar  transducers  is  critical.  Due  to  the  fact  that  each  transducer 
radiates  in  close  proximity  to  one  another,  mutual  interference  occurs.  The  technique 
employed  in  common  practice  to  model  this  array  element  interaction  [9,  Chp  8]  is  not 
applicable  when  the  array  elements  are  closely-packed.  This  led  to  the  emergence  of  a 
research  program  to  develop  the  so-called  Transition,  or  T-matrix  method  for  modeling 
the  radiation  from  an  arbitrarily  densely-packed  array  of  sonar  transducers  [1,3]. 

The  T-matrix  relates  the  incident  pressure  applied  on  a  particular  transducer  to  its 
scattered  pressure  at  a  specific  harmonic  component.  The  scattering  properties  can  be 
determined  through  Finite  Element  Modeling  (FEM),  where  a  known  pressure  field  is 
imposed  on  the  transducer,  in  this  case,  represented  by  a  thin  spherical  shell.  The  T- 
matrix  of  a  single  transducer  will  subsequently  be  incorporated  into  the  calculation  of  the 
scattered  pressures  from  an  array  of  sonar  transducers. 

The  analytical  T-matrix  can  be  determined  from  the  modal  mechanical  impedance 
of  a  spherical  shell,  which  relates  the  pressures  and  the  displacement  components.  These 
analytical  values  form  the  basis  for  validation  of  the  modeled  results.  In  a  previous 
attempt  using  the  ATILA  finite-element  code,  the  computed  T-matrix  was  deemed 
unsatisfactory  in  some  aspects  [2].  The  discrepancies  could  be  due  to  the  ineffectiveness 
of  the  radiation  boundary  condition  and  the  coarse  mesh  built  for  the  model. 

B.  PURPOSE 

This  thesis  investigates  the  use  of  COMSOL  Multiphysics  finite-element  code  for 
the  computation  of  the  T-matrix  for  a  thin  spherical  shell.  A  3D  model  of  a  water-loaded 
spherical  thin  steel  shell  with  an  outer  radius  of  0.5m  and  a  thickness  of  1cm  was 
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constructed,  employing  several  fine  meshes.  The  scattering  T-matrix  elements  for  this 
shell  were  computed  at  a  frequency  of  474  Hz  (i.e.,  ka=  1).  Two  acoustic-structure 
interaction  types  of  models  were  created,  one  employing  2-D  thin  shell  elements  to  model 
the  physical  shell  and  one  in  which  3-D  elastic  elements  were  used  to  model  the  physical 
shell. 

In  the  model  calculations,  an  incident  wave  distributed  in  angle  as  one  spherical 
hannonic  component  was  applied  and  the  corresponding  scattered  waves,  in  tenns  of  the 
spherical  harmonic  components,  were  evaluated.  Scattering  coefficients  for  the  zero 
(monopole),  first  (dipole)  and  second  (quadrupole)  order  modes  were  determined.  The 
computed  T-matrix  was  validated  with  reference  to  analytical  values. 

C.  THESIS  OUTLINE 

The  subsequent  chapters  of  the  thesis  are  outlined  as  follows: 

•  Chapter  II  covers  two  parts  of  the  theory.  The  first  part  entails  the  T- 
matrix  theory  involving  the  necessary  equations  that  lead  to  the 
formulation  for  a  single  transducer  and  an  array  of  transducers.  The 
second  part  illustrates  the  theory  of  forced  vibration  of  a  spherical  shell 
and  a  series  of  manipulations  used  to  derive  the  analytical  T-matrix. 

•  Chapter  III  introduces  the  features  in  the  COMSOL  Multiphysics  Finite- 
Element  code  and  presents  the  3D  model  of  the  water-loaded  thin 
spherical  shell. 

•  Chapter  IV  presents  the  T-matrix  results  detennined  by  COMSOL 
modeling  and  the  percentage  of  error  with  respect  to  the  analytical  values. 

•  Chapter  V  presents  the  conclusions  of  the  thesis  and  suggestions  for  future 
work. 
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II.  THEORY 


A.  THE  “T-MATRIX”  APPROACH 

The  “T-matrix”  relates  the  incident  and  scattered  waves  from  an  obstacle  or 
scatterer.  From  the  free- Held  radiation  and  (short-circuit)  scattering  properties  of  a  single 
active  sonar  transducer,  as  represented  in  the  T-matrix,  one  can  compute  the  radiation 
properties  of  an  array  of  arbitrarily  densely  and  randomly  packed  sonar  transducers  [1]. 
The  relationship  between  the  incident  pressure,  P1 ,  and  the  scattered  pressure,  Ps ,  on  a 
single  element  is  represented  by: 

{p'HTKpi}  (1) 

where  {  }  represents  a  column  of  coefficients  and  [T]  is  the  T-matrix. 

1.  Single  Transducer 

For  a  particular  transducer  submerged  in  fluid,  an  incident  wave  is  represented  by 
a  standing  spherical  wave.  It  is  comprised  of  a  discrete  basis  set  of  (unnormalized) 
spherical  hannonics,  with  an  e“ot  harmonic  time  dependence  assumed,  and  is  expressed 
in  the  fonn  of: 

pi  M.40 = S"  »SL,p-j.(kr>p"<cos0)ei”*  o) 

The  scattered  pressure  is  similarly  represented  by: 

P‘  (r.  0, <l>) =  Z"  0  ZL„  O™  (kr)p;  (cos  6)e™*  (3) 

where: 

n,m,9,(j)  :  spherical  harmonics  indices  of  arguments 

P|m  :  coefficient  of  the  incident  pressure 

jn(kr)  :  nth  order  spherical  Bessel  function  of  the  first  kind 
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Pnm  (cos  0)  :  associated  Legendre  function  of  the  first  kind 

h(n2)  :  nth  order  spherical  Hankel  function  of  the  second  kind 

2.  Array  of  Transducers 

In  an  array  of  active  transducers,  where  their  mutual  field  interactions  are  to  be 
taken  into  account,  the  T-matrix  formulation  will  have  to  involve  the  use  of  the  spherical 
addition  theorem.  Considering  an  array  of  two  transducers,  illustrated  in  Figure  1  [2,  Chp 
2],  whereby  their  radiated  fields  interact  with  each  other,  the  formulation  is  as  follows: 


TjG12(A2  +B2)  -  Aj 

(4) 

T2G21(Aj  +  B,)  =  a2 

(5) 

where: 

Gi,  : 

translation  formula  from  system  j  to  system  i 

(by  the  spherical  addition  formula) 

Tk  : 

T-matrix  for  the  k1h  transducer  modeled 

Bj  and  B2  : 

radiated  amplitudes  of  the  two  transducers 

Aj and  A2  : 

scattering  pressure  amplitudes  of  the  two  transducers 

4 


Sphere  1 


Sphere  2 


zafleidpolnl 


Figure  1.  Interaction  between  an  array  of  two  transducer  (From  Ruiz,  1994)  [2] 


Putting  into  matrix  form: 


'  I  -T.G.P 

t,g12b2s 

V-T2G2i  I  y 

V^2y 

J2G21Bly 

After  simplification,  the  resulting  equation  becomes: 

(1-3) 


PC 

=  d 

V^-2  ) 

Ik 

where: 


'  0 

VT2G21 


t,g12 

0 


x 

) 


(6) 


(7) 


(8) 


Upon  the  detennination  of  Al  and  A,  ,  the  total  pressure  at  any  field  point  can  be 
obtained  by  summation  over  the  entire  array  element: 

PC. ».  ♦>  =  ZZZ  ZL.  [A-  +  Bi™  »j)  C(cos9)e™»  (9) 

j 
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As  an  example,  this  methodology  enables  the  determination  of  the  source  levels 
for  an  array  of  16  closely  packed  spherical  shell  transducers  operating  at  474Hz  (ka=  1 )  as 
shown  in  Figure  2  [1]: 


Orientation  of  the  spherical  shells 
in  the  16  element  array 


Figure  2.  Array  of  16  closely  packed  spherical  shell  transducer  (From  Baker,  1998) 

[1] 


B.  FORCED  VIBRATION  OF  A  THIN  SPHERICAL  SHELL 

We  follow  the  analysis  by  Junger  and  Feit  of  the  forced  vibration  of  a  thin 
spherical  shell  [4,  Chp  7],  In  the  study  of  acoustic-structure  interaction  problem,  they 
adopted  Loves  ’s  assumptions  for  the  equations  of  motion  of  thin  elastic  shells: 

•  The  thickness  of  a  shell  is  small  compared  with  the  smallest  radius  of  the 
curvature  of  the  shell. 

•  The  displacement  is  small  in  comparison  with  the  shell  thickness. 

•  The  transverse  nonnal  stress  acting  on  planes  parallel  to  the  shell  middle 
surface  is  negligible. 

•  Fibers  of  the  shell  normal  to  the  middle  surface  remain  so  after 
deformation  and  are  themselves  not  subject  to  elongation. 

In  addition,  the  shell  is  thin  enough  and  the  frequency  low  enough  to  ignore  flexural 
stresses  as  compared  with  membrane  stresses  [4,  Chp  11]. 
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As  a  result,  the  displacements  of  the  spherical  shell  can  be  perceived  as  a 
function  of  the  shell’s  middle  surface  deflection,  shown  in  Figure  3  [5], 


In  view  of  nontorsional  axisymmetric  motions  of  a  spherical  shell,  the  tangential 
component  vanishes  in  (j> ,  the  azimuthal  direction.  The  existing  components  are  w,  the 
mid-surface  radial  motion  and  u,  the  tangential  motion  in  the  direction  of  increasing  0 , 
and  both  are  independent  of  (j) .  A  further  optimization  [2,  Chp  2]  reduced  the  equation  to 
contain  only  the  w  component.  The  resulting  equation  of  interest  is,  therefore,  the  modal 
mechanical  impedance,  given  by  [4,  Chp  7]: 

Z  n=T^SS-  (10) 
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The  respective  incident  pressure,  scattered  pressure  and  radial  displacement  can  be 
expressed  as: 


(12) 

Z,„R..hi2|(kr)p:(cose)e”* 

(13) 

Z„n.w.n.C(cose)e"» 

(14) 

where  Inm  and  Rnm  are  the  incident  and  reflected  modal  amplitudes  respectively,  and 
wnm  is  the  modal  normal  displacement  of  the  shell.  The  pressure  at  the  shell  surface  is: 

P  =  P‘+PS  (15) 


Expressing  the  modal  mechanical  impedance  as  [4,  Chp  7]: 


Z 


n 


hCpPs 

r  J 

[q2  -2(1  +  u)](q2  +l-u-Zn)-Zn  (l  +  n)2 

a Q 

Q2  +l-u-A,n 

(16) 


where 


K  =n(n  +  l) 


C 


co  :  angular  frequency 

h  :  shell  thickness 

Cp  :  plate  phase  velocity 

o  :  Poisson’s  ratio 

E  :  Young’s  modulus 

ps  :  shell  density 

a  :  shell  radius 

Arranging  Equations  (10)  and  (11)  yields: 
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(17) 


T„J„  (ka)  +  Rnmhi2)  (ka)  =  -icoZnw 


Another  relationship  between  Inmj  Rnmand  wnm  can  be  established  using  Euler’s 


equation. 


<9p .  _  S2w 

~dr'I=a~~Pf~dtr 


(18) 


where  pf  is  the  density  of  the  fluid,  and  with  the  assumption  of  no  cavitation  on  the 
surface  of  the  shell, 

I„mjn(ka)  +  Rnmhi2)'(ka)  =  copfCfwnm  (19) 


Eliminate  wnm  by  equating  Equations  (17)  and  (19),  the  relationship  between  Inm  and 
Rnm  is  simplified  into: 


iz„j,(ka>  +  PfC,j„(ka)  1 

[iZnhf(ka)  +  prCfh®(ka)J 


Equation  (20)  will  be  used  to  derive  the  analytical  results  of  the  T-matrix  in  which  each 
value  of  Rnm  corresponds  to  the  diagonal  element  of  the  matrix  when  Ln  is  set  to  one. 


The  T-matrix  format  for  the  hannonics  up  to  the  second  order  (n=2)  is  as  follows: 


T 


R00  0  0  0  0  0  000 

0  R,_,  0  0  0  0  0  0  0 

0  0  R10  0  0  0  000 

0  0  0  Rn  0  0  000 

0  0  0  0  R2_2  0  0  0  0 

0  0  0  0  0  R2_!  0  0  0 

ooooo  or20oo 

00000  0  0  R210 

OOOOO  0  00  R22 


s-' 


Zero  Order  First  Order  Second  Order 

(Monopole)  (Dipole)  (Quadrupole) 


(21) 
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Due  to  spherical  symmetric,  the  scattered  wave  should  contain  the  same  harmonic 
components  as  the  incident  wave.  In  other  words,  the  resulting  T-matrix  should  be 
diagonal  and  the  non-diagonal  values  should  be  zero. 
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III.  COMSOL  MULTIPHYSICS  FINITE-ELEMENT  CODE 


This  chapter  introduces  the  useful  features  of  the  Acoustic  Module  in  COMSOL 
Multiphysics  [6]  and  describes  the  design  of  the  3D  spherical  thin  shell  model  and  its 
essential  governing  equations.  For  the  purpose  of  achieving  accurate  results,  the  model 
simulation  was  attempted  using  two  different  “physics  modes”  -  the  Acoustic-Shell 
Interaction  (ACSH)  and  the  Acoustic-Structure  Interaction  (ACSI)  modes.  The  key 
difference  between  these  two  physics  modes  is  the  theory  that  is  governing  the  solid.  The 
ACSH  is  governed  by  thin  shell  approximation  whereas  the  ACSI  employs  the  full  3D 
elastic  theory.  The  rest  of  the  chapter  will  cover  the  fonnulation  of  equations  for  the 
incident  standing  wave  as  well  as  the  computation  of  the  T-matrix  elements  with  the 
resulting  scattered  pressure  obtained  in  COMSOL. 

A.  INTRODUCTION  TO  COMSOL  ACOUSTIC  MODULE 

The  COMSOL  Multiphysics  Finite-Element  Code  is  capable  of  solving  real-time 
physics  problem  with  precision.  It  contains  a  wide  variety  of  modules  for  different 
“physics  modes”  that  enables  the  user  to  define  the  conditions,  physical  effects  and 
mathematical  equations  that  are  essential  in  the  modeling.  The  Acoustic  Module  featured 
in  COMSOL  Multiphysics  Version  4.2a  contains  a  collection  of  physics  interfaces  that 
are  adapted  to  a  broad  category  of  acoustics  simulations  in  fluids  and  solids.  It  supports 
transient,  eigenfrequency,  frequency  domain,  mode  analysis,  and  boundary  mode  analysis 
in  pressure  acoustics  in  compressible  and  irrotational  velocity  fields  [6].  In  this  thesis,  the 
two  physics  modes  used  are  as  follows: 

1.  Acoustic-Structure  Interaction  Physics  Mode 

This  is  a  multiphysics  mode  where  the  fluid’s  pressure  causes  a  fluid  load  on  the 
solid  domain,  and  the  structural  acceleration  affects  the  fluid  domain  as  a  normal 
acceleration  across  the  fluid-solid  boundary.  In  the  frequency  domain,  the  harmonic 
sound  waves  in  the  fluid  domain  are  governed  by  the  Helmholtz  equation: 
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V- 


(22) 


1  /  ^  co2p 

- (Vp  +  q) - 4  =  0 

l  Pc  )  PcC“ 

where  the  acoustic  pressure  is  p  =  p0elwt  in  (N/m2)  and 
pc  :  fluid  density  (kg/m3) 

q  :  optional  dipole  source  (m/s2) 

co  :  angular  frequency  (rad/s) 

c  :  speed  of  sound  (m/s) 

The  default  boundary  conditions  that  influence  the  fluid-solid  interaction  are  the 
pressure  load  (force/unit  area)  by  the  fluid  acting  on  the  solid, 

F  =  -nsp  (23) 

where  nsis  the  outward-pointing  unit  normal  vector  seen  from  inside  the  solid  domain, 
and  the  normal  structural  acceleration  at  the  boundary  [6], 

a„=ns-utt  (24) 

where  utt  is  the  second  derivative  of  the  structural  displacement  with  respect  to  time. 

2.  Acoustic-Shell  Interaction  Physics  Mode 

This  mode  uses  the  frequency  domain  interface  feature  and  a  thin  shell 
approximation  to  model  the  interactions  between  a  deformable  shell  and  the  acoustic 
pressure  in  a  fluid  domain.  The  shell  thickness  is  defined  by  the  user.  The  boundary 
condition  to  model  any  deformable  shell  boundary  can  be  chosen  to  be  an  exterior  shell 
or  an  interior  shell.  For  the  exterior  shell,  only  one  side  of  the  boundary  is  adjacent  to  the 
fluid  domain  and  the  normal  acceleration  is  given  by: 

(  i 

-v - (vP-q)  =ns'Utt  (25) 

l  Po  J 

and  the  pressure  load  on  the  shell  is  given  by  Equation  (23). 
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As  for  the  interior  shell,  both  sides  of  the  boundary  interact  with  the  fluid  domain. 
The  normal  acceleration  on  each  side  is  akin  to  Equation  (25)  and  the  pressure  load  on 
the  shell  is  [6], 

Fp  =-(nsp)i -(nsp)2  (26) 

B.  THE  3D  MODEL 

A  three  dimensional  model  of  a  thin  spherical  steel  shell  was  created  in  COMSOL 
as  shown  in  Figure  4.  The  dimension  of  the  thin  shell  is  based  on  the  size  of  a  typical  low 
frequency  transducer  at  ka~  1.  As  mentioned  by  Timme,  Young  and  Blue  [9,  Chp  1],  for 
the  case  of  ka  <  1,  the  transducer’s  radiated  acoustic  power  is  deemed  inefficient.  In  the 
present  model,  the  approach  to  obtain  ka=  1  is  to  operate  the  transducer  at  474Hz,  which 
is  close  to  the  natural  frequency  of  submerged  spherical  shell  suggested  by  Junger  and 
Feit  [4,  Chp  9]. 


Figure  4.  3D  Model  of  a  thin  spherical  shell  submerged  in  fluid 


1. 


Model  Parameters 


The  material  of  the  thin  shell  is  steel  and  the  parameters  employed  in  the  model 
are  defined  as  follows: 

a  Shell  outer  radius  0.5m 
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h  Shell  thickness  0.01m 

R  Outer  radiation  boundary  of  the  external  fluid  2.5m 

E  Young’s  modulus  (steel)  215GPa 

v  Poisson’s  ratio  (steel)  0.33 

/  Operating  frequency  474Hz  (ka=l) 

ps  Shell  density  (steel)  7500kg/m3 

pf  Fluid  density  1000kg/m3 

Cf  Sound  speed  in  fluid  1490m/s 


2.  Radiation  boundary 

A  spherical  wave  radiation  boundary  condition  is  specified  over  the  fluid 
outer  boundary  to  absorb  any  outgoing  waves.  Alternatively,  a  perfectly  matched  layer 
(PML)  can  be  used  to  exhibit  the  same  feature. 

3.  Meshing 

Tetrahedral  elements  were  employed  to  create  the  mesh.  According  to 
COMSOL  guideline  [6],  the  mesh  element  should  be  no  greater  than  one-fifth  of  the 
operating  wavelength.  For  the  operating  frequency  of  474Hz  employed,  the  wavelength 

in  fluid  is  3.14m  and  the  maximum  element  size  should  be  —  A  =  0.63m  or  smaller.  The 

5 

meshed  model  shown  in  Figure  5  has  a  so-called  “finer”  (COMSOL  tenninology)  mesh 
for  the  fluid  domain  and  an  “extremely  fine”  mesh  for  the  thin  shell. 
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Fluid  domain: 

Max  element  size:  0. 175m  (0.056  X ) 
Min  element  size:  0.0075m  (0.0024  X ) 

Thin  shell: 

Max  element  size: 

0.1m  (10  x  thickness) 

Min  element  size: 

0.001m  (0.1  x  thickness) 


Figure  5.  Model  with  finer  mesh  for  the  fluid  domain  and  extremely  fine  mesh  for  the 

thin  shell 

4.  Incident  wave  onto  shell  surface 

The  desired  effect  is  to  impose  a  spherical  standing  wave  acting  upon  the 
shell  surface  as  a  force  per  unit  area.  The  equation  for  the  incident  pressure  was 
expressed  in  Equation  (2)  and  is  reproduced  as  Equation  (27).  As  the  “Incident  pressure” 
option  in  COMSOL  only  acts  upon  the  radiation  boundary,  the  decision  was  to  use  the 
“Background  pressure”  option  to  specify  the  spherical  standing  wave  throughout  the 
volume  of  the  fluid,  to  achieve  the  desired  effect.  The  equation  of  the  incident  pressure 
will  be  elaborated  in  the  next  section. 

C.  COMSOL  MODELING  FOR  T-MATRIX  DETERMINATION 

The  two  important  processes  for  the  computation  of  the  T-matrix  are  specifying 
the  incident  standing  wave  (spherical  Bessel  functions)  and  the  post  evaluation  of  the 
scattered  pressure  (Hankel  functions)  over  the  shell  surface  to  obtain  the  elements  of  the 
T-matrix.  The  workflow  of  the  thin  shell  modeling  in  COMSOL  is  shown  in  Figure  6. 
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Incident  standing  wave 
(Spherical  Bessel  Function) 


Outgoing  wave 
(Hankel  Function) 


Post 

Evaluation 


Thin  Shell 


COMSOL  Modeled  Interaction:  ACSI/ACSH 


T-MATRIX 


Figure  6.  Workflow  of  COMSOL  Modeling 

Recall  that  the  incident  and  scattered  waves  are  taken  to  be  Equations  (2)  and  (3), 
required  here: 

P'(r.^)  =  Z:»ZL„Ci(kr)C(cose)e-  (27) 


P‘  (> r.  < 9. 40  =  I"  „  2L_„ ' Plh®  (kr)p,r  (cos  0)e" 


Accounting  for  scattering  up  through  quadrupolar  (n=0,  1,  2;  m=-n,  n),  the  T- 
matrix  between  the  incident  and  scattered  waves  (Equation  (1))  takes  the  following  form: 


Each  column  of  the  T-matrix  is  obtained  by  setting  each  Rm  equal  to  one,  one  at  a  time, 
while  setting  all  the  other  P'rn  to  zeros,  and  employing  orthogonality  of  the  spherical 
harmonics  to  extract  the  set  of  Pnsm  for  all  n  and  m  of  interest  (n=0,  1,  2).  For  instance,  by 
setting  only  Pq0  =  1  and  the  others  to  zero,  the  first  column  of  the  T-matrix  will  be 
equaled  to  the  corresponding  Pnsm  column.  The  procedure  is  as  follows: 
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Note,  the  product  Pnm(cos0)eim<l>  in  p1  and  ps  is  related  to  the  usual  (normalized) 

spherical  harmonic  representation  of  angular  dependence  in  spherical  coordinates  by  [c.f., 
7,  Chp  18] 

Y”(0,+)  =  C„P;  (cos0)exp(im+) ,  (30) 


where 


C„ 


!-l)m 

2n  +  l  (n-m)! 

1/2' 

(  V 

An  (n  +  m)! 

(31) 


Normalized  in  this  way,  the  spherical  harmonics  have  the  orthonormality  property: 

1  2  7T  * 

J  J[Y.”(e,+)]  Y„"'(e,*)d*d(cos0)  =  81„.5_..  (32) 

-1  0 


In  terms  of  the  Y”  (0,  <j)) ,  any  function  of  angle  can  be  expanded  as  [7,  Chp  18]: 

f(e.*)=ZLIL,^Y;(e,*),  (33) 

where 

1  2  71  * 

a-  =  ]  |[y;(0,P]  f(e,«|))dit.d(cos0)  =  CmIm,  (34) 

-1  0 


1  271 

1  nm  =  |  J  P™  (cos  9)  exp  (-im(|))f  (0, 4>)  dc))d(cos  0) .  (35) 

-1  0 


Comparing  Equations  (28)  and  (33),  setting 

Z:=.ZL„4,«Y;(0,4,)  =  XZZL„Cl>!,2’(kr)P;(cos0)e' 


im<j) 


(36) 


and  eliminating  the  common  terms,  the  desired  expression  for  scattered  pressure 
amplitude  Pnsm  is, 
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(37) 


PI  = 


nm  nm 


Cl  •  I„ 


nm  nm 


h«(kr)  h«'(kr)  ' 


Equation  (37)  will  be  used  in  the  post  evaluation  of  the  modeled  result  Inm.  The  required 

integral,  is  computed  on  the  shell  surface,  within  COMSOL  Multiphysics,  using  the  nodal 
values  of  the  scattered  pressure. 

Pnsmis  equivalent  to  the  T-matrix  element  Rnm  of  spherical  harmonics  n  and  m.  Since  the 
incident  pressure  amplitudes  are  always  set  equal  to  one,  in  only  one  spherical  harmonic 
at  a  time,  the  values  of  Pi  so  determined  are  the  T-matrix  element  columnar  values. 
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IV.  RESULTS 


This  chapter  presents  the  T-matrix  elements  computed  for  the  specified  spherical 
thin  shell,  obtained  using  COMSOL  Multiphysics.  We  consider  our  goal  as  to  achieve  an 
accuracy  of  0.1%  for  the  scattering  properties  of  a  single  sonar  transducer.  We  base  this 
goal  on  the  supposition  that,  in  a  worst  case,  there  might  occur  as  much  as  90% 
destructive  interference  between  elements  in  an  array.  Then,  if  we  desire  a  final  array 
model  accuracy  of  1%  (i.e.,  O.ldB),  we  need  an  accuracy  of  0.1%  for  the  computation  of 
single-element  scattering. 

For  a  start,  we  compute  the  T-matrix  of  a  rigid  sphere  in  fluid,  a  less  complex 
model,  to  examine  and  become  familiar  with  the  functionality  of  COMSOL.  This  model 
was  also  used  to  detennine  the  best  fluid  mesh  parameters.  We  then  consider  the  analysis 
in  using  thin  shell  elements  to  model  the  steel  shell,  to  examine  the  ability  and  accuracy 
of  the  Acoustic-shell  Interaction  physics  mode  to  model  single-element  scattering. 
Finally,  we  consider  the  analysis  in  using  true  3D  elastic  finite  elements  to  model  the 
shell  to  validate  the  ability  of  the  Acoustic-structure  Interaction  physics  mode  to  model 
the  scattering.  The  T-matrix  obtained  from  each  simulation  will  be  verified  against  an 
analytical  computation  of  the  T-matrix. 

In  all  subsequent  calculations  of  the  T-matrix  elements,  we  restrict  ourselves  to 
monopolar  (n=0),  dipolar  (n=l)  and  quadrupolar  (n=2)  terms.  (See  Equation  (21)) 

A.  T-MATRIX  OF  A  RIGID  SPHERE  IN  FLUID 

The  purpose  of  this  section  is  to  examine  the  functionality  of  COMSOL  by 
modeling  a  fluid  loaded  rigid  sphere.  This  model  was  simulated  in  the  COMSOL 
Multiphysics  Acoustic-shell  Interaction  (ACSH)  physics  mode  with  the  geometry  and 
physical  properties  listed  in  Chapter  III,  Section  B,  except  that  the  material’s  density  and 
Young’s  modulus  were  increased  to  106  times  the  original  properties,  to  simulate  a  rigid, 
stationary  sphere. 
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As  illustrated  in  Chapter  III,  Section  C,  an  incident  standing  wave  (spherical 
Bessel  function)  as  of  Equation  (27)  is  imposed  onto  the  object.  The  scattered  pressure 
(Hankel  function)  resulted  from  the  interaction  is  represented  by  Equation  (28).  As  such, 
the  T-matrix  diagonal  elements  can  be  determined  by  evaluating  the  scattered  pressure 
over  the  sphere  surface  using  Equation  (37). 

The  analytical  T-matrix  values  for  the  rigid  sphere  are  presented  in  Table  1.  These 
were  computed  using  Equation  (38),  which  was  evaluated  from  Equation  (20)  by  taking 
the  limit  Z  — »  oo  : 

n 

R  _  Jnjka)  _  (n  + 1)  jn+1  (ka)-n  •  jn_,  (ka) 

hf(ka)-n.h^(ka)-(n  +  l)h<2>(ka)  1  > 

Table  2  presents  two  sets  of  model  results,  the  first  obtained  using  a  normal  mesh 
for  both  the  fluid  and  the  solid  surface,  and  the  second  obtained  using  more  refined 
meshes  for  the  fluid  and  solid  surface.  It  only  presents  the  diagonal  T-matrix  elements; 
the  non-diagonal  elements  should  equal  zero  and  will  be  considered  in  a  later  section. 


^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Roo 

-4.5351E-02 

2.0807E-01 

2.1296E-01 

1.0230E+02 

Rio 

-1.1437E-02 

-1.0633E-01 

1.0694E-01 

-9.6139E+01 

R20 

-1.4876E-04 

-1.2196E-02 

1.2197E-02 

-9.0699E+01 

Table  1.  Analytical  T-matrix  for  a  rigid  sphere 
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Normal  mesh  for  fluid  and  solid 

^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

Roo 

-4.5050E-02 

2.0767E-01 

2.1250E-01 

1.0224E+02 

-0.215 

-0.0564 

Rio 

-1.1400E-02 

-1.0597E-01 

1.0658E-01 

-9.6140E+01 

-0.339 

-0.0009 

R20 

-1.4573E-04 

-1.2070E-02 

1.2071E-02 

-9.0692E+01 

-1.032 

0.0070 

Refined  Mesh  (Fluid:  extra  fine  meshes  and  solid:  extremely  fine  meshes) 

^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

Roo 

-4.5215E-02 

2.0779E-01 

2.1265E-01 

1.0228E+02 

-0.142 

-0.0199 

Rio 

-1.1413E-02 

-1.0620E-01 

1.0681E-01 

-9.6134E+01 

-0.123 

0.0053 

R20 

-1.4686E-04 

-1.2155E-02 

1.2156E-02 

-9.0692E+01 

-0.335 

0.0066 

Table  2.  COMSOL:  T-matrix  for  a  rigid  sphere 


The  T-matrix  elements  for  the  rigid  sphere  were  evaluated  on  the  sphere  surface. 
For  the  nonnal  mesh,  the  element  size  is: 

Fluid: 

Maximum  element  size  =  0.5  m  (0.16X) 

Sphere  surface: 

Maximum  element  size  =  0.2  m  (0.064^) 

As  for  the  refined  mesh,  the  element  size  is: 

Fluid  (extra  fine  mesh): 

Maximum  element  size  =  0.175  m  (0.056A.) 

Sphere  surface  (extremely  fine  mesh): 

Maximum  element  size  =  0.1  m  (  0.032a) 

From  Table  2,  quadrupole  R20  is  observed  to  be  the  worst  case,  with  the  highest 

magnitude  error.  After  a  refined  mesh  was  introduced,  its  magnitude  error  with  respect  to 

the  analytical  values  shows  an  improvement  from  -1.032%  to  -0.335%.  In  view  of  the 

goal  of  0.1%  accuracy,  these  results  are  considered  marginally  satisfactory. 
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B.  T-MATRIX  OF  A  THIN  SPHERICAL  SHELL  USING  THIN  SHELL 
APPROXIMATION 

1.  Analytical  T-matrix 

The  analytical  values  are  determined  by  Equation  (20)  using  thin  shell 
approximation.  The  T-matrix  diagonal  elements  are  shown  in  Table  3. 


J^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Roo 

Rio 

R20 

-1.3696E-02 

-6.0902E-03 

-6.2949E-04 

1.1623E-01 

7.7802E-02 

-2.5082E-02 

1 . 1 703E-01 

7.8039E-02 

2.5090E-02 

9.6721E+01 

9.4476E+01 

-9.1438E+01 

Table  3.  Analytical  T-matrix  for  a  thin  shell  at  n=0,  1  and  2 


2.  COMSOL  Results  for  Thin  Shell  Approximation  Modeling 
a.  Diagonal  elements  of  T-matrix 

The  thin  spherical  shell  was  modeled  in  the  COMSOL  Multiphysics 
Acoustic-shell  Interaction  (ACSH)  physics  mode  with  extra  fine  meshes  for  the  fluid 
domain  and  extremely  fine  meshes  for  the  shell  surface.  In  this  mode,  the  shell  is 
modeled  using  thin  shell  finite  elements,  that  is,  the  shell  is  modeled  with  a  thin  shell 
approximation  [4,  Chp  7].  The  mesh  element  sizes  are  similar  to  the  refined  mesh  of  the 
rigid  sphere  in  previous  section. 

The  mesh  of  the  model  consisted  of  214299  elements  and  298606  degrees 
of  freedom  (DOF).  It  was  computed  over  an  average  time  frame  of  12  minutes.  The 
diagonal  T-matrix  elements  for  this  model,  evaluated  on  the  shell  surface,  are  presented 
in  Table  4.  The  last  two  columns  show  the  magnitude  and  phase  errors  relative  to  the 
analytical  T-matrix. 
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^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

Roo 

-1.3693E-02 

1.1623E-01 

1 . 1 703E-01 

9.6719E+01 

0.0010 

-0.0019 

Ri,-i 

-6.0884E-03 

7.7799E-02 

7.8037E-02 

9.4475E+01 

-0.0027 

-0.0012 

Rio 

-6.0831E-03 

7.7767E-02 

7.8004E-02 

9.4473E+01 

-0.0449 

-0.0032 

Rn 

-6.0884E-03 

7.7799E-02 

7.8037E-02 

9.4475E+01 

-0.0027 

-0.0012 

R2,-2 

-6.2334E-04 

-2.501  IE-02 

2.5019E-02 

-9.1428E+01 

-0.2818 

0.0100 

R2,-i 

-6.2475E-04 

-2.5039E-02 

2.5047E-02 

-9.1429E+01 

-0.1706 

0.0084 

R20 

-6.2540E-04 

-2.5053E-02 

2.5061E-02 

-9.1430E+01 

-0.1133 

0.0077 

R2I 

-6.2475E-04 

-2.5039E-02 

2.5047E-02 

-9.1429E+01 

-0.1706 

0.0084 

R22 

-6.2335E-04 

-2.501  IE-02 

2.5019E-02 

-9.1428E+01 

-0.2818 

0.0100 

Table  4.  COMSOL:  Diagonal  T-matrix  elements  for  a  thin  shell  (fluid  domain:  extra  fine 

meshes  and  shell  surface:  extremely  fine  meshes) 


From  the  results  presented  in  Table  4,  three  observations  were  made.  Firstly,  the 
magnitude  error  for  R0o,  Rio  and  R20,  ascertained  that  the  results  are  in  better  agreement 
than  those  of  the  rigid  sphere.  The  discrepancy  is  unknown  at  this  juncture  and  remains  to 
be  investigated. 

The  second  observation  deduced  that  the  results  are  almost  independent  of  the 
azimuthal  mode  number,  m,  as  they  should  be  for  the  case  of  an  axisymmetric  scatterer. 
However,  for  the  modeling  in  COMSOL,  the  axisymmetric  characteristics  could  have 
been  slightly  affected  by  the  differing  mesh  sizes  built  on  the  model. 

The  third  observation  arises  from  the  comparison  with  the  T-matrix  determined 
by  ATILA  Finite-Element  Modeling,  with  its  results  presented  in  Table  5  [2,  Chp  3],  The 
model  was  composed  of  2746  nodes  and  288  mesh  elements.  The  number  of  mesh 
elements  is  nearly  a  thousand  times  lesser  than  the  COMSOL  mesh  elements  of  214299. 
It  can  be  seen  from  the  last  two  columns  that  a  high  magnitude  error  of  up  to  -38% 
occurred  in  the  quadrupole  scattering  computation.  It  is,  therefore,  evident  that  a  mesh 
refinement  in  COMSOL  modeling  has  yielded  more  accurate  results. 
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I^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

o 

o 

04 

-1.2624E-02 

1. 1429E-01 

1.1499E-01 

9.6303E+01 

-0.91 

-0.36 

Ri,-i 

-6.8902E-03 

8.2376E-02 

8.2664E-02 

9.4781E+01 

6.49 

0.33 

Rio 

-5.7422E-03 

8.3113E-02 

8.3311E-02 

9.3952E+01 

7.33 

-0.50 

Rn 

-6.8668E-03 

8.2363E-02 

8.2649E-02 

9.4766E+01 

6.47 

0.31 

R2.-2 

-2.8487E-03 

-1.5203E-02 

1.5468E-02 

-1.0061E+02 

-37.8 

-9.19 

R2.-1 

-2.7899E-03 

-1.5159E-02 

1.5414E-02 

-1.0043E+02 

-38.0 

-9.00 

R20 

-2.4889E-03 

-1.5287E-02 

1.5488E-02 

-9.9247E+01 

-37.7 

-7.82 

R21 

-2.7847E-03 

-1.5164E-02 

1.5418E-02 

-1.0041E+02 

-38.0 

-8.98 

R22 

-2.8227E-03 

-1.5182E-02 

1.5442E-02 

-1.0053E+02 

-37.9 

-9.11 

Table  5.  ATILA:  Diagonal  T-matrix  elements  for  a  thin  shell  [2,  Chp  3] 


b.  Non-diagonal  elements  of  T-matrix 

The  non-diagonal  element  is  a  measure  of  the  T-matrix  efficacy.  A 
significant  value  indicates  an  existence  of  a  numerical  “leakage”  in  the  model 
computation.  In  the  T-matrix  detennined  by  ATILA,  there  are  two  significant  leakages 
[2,  Chp  4]. 

The  T-matrix  in  Table  6  presents  the  non-diagonal  elements  relative  to  the 
diagonal  elements  in  the  same  column.  This  normalization  was  chosen  because  the  T- 
matrix  elements  are  computed  one  column  at  a  time.  For  instance,  the  non-diagonal  value 
for  Rio  (third  column)  relative  to  R0o  (first  row)  corresponds  to  4.98E-5. 

The  largest  three  non-diagonal  values  are  1.63E-03  (monopole  to  dipole), 
1.63E-03  (dipole  to  quadrupole)  and  1.16E-02  (quadrupole  to  quadrupole),  as  highlighted 
in  Table  6.  However,  these  values  are  not  considered  to  be  significantly  causing  a  leakage 
in  the  model.  As  a  comparison  to  ATILA  modeling,  the  reduction  in  leakage  could  have 
been  due  to  the  finer  mesh  employed  in  the  COMSOL  modeling. 
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Scattered  Pressure 


Incident  Pressure 


Roo 

R1-1 

Rio 

Rn 

R2,-2 

R2,-l 

R20 

R21 

R22 

Roo 

~  1 

1.80E-04 

4.98E-05 

3.61E-04 

4.76E-07 

6.34E-06 

5.78E-05 

4.24E-05 

1.32E-05 

Ri,-i 

1.63E-03 

1 

6.16E-05 

2.40E-05 

6.77E-07 

8.13E-06 

3.14E-04 

2.22E-05 

9.27E-06 

Rio 

2.41E-04 

3.07E-05 

1 

6.19E-05 

1.61E-06 

1.40E-05 

5.98E-05 

8.05E-05 

3.91E-05 

Rn 

8.15E-04 

6.04E-06 

3.05E-05 

1 

1.51E-07 

1.83E-06 

1.56E-04 

2.42E-05 

1.10E-05 

R2,-2 

2.79E-04 

1.34E-04 

1.80E-04 

6.88E-05 

1 

3.67E-05 

8.61E-05 

3.70E-04 

1.16E-02 

R2,-l 

8.54E-04 

1.34E-04 

4.48E-04 

1.60E-04 

9.16E-06 

1 

2.08E-04 

4.62E-04 

3.70E-04 

R20 

7.85E-04 

8.13E-04 

2.19E-04 

1.63E-03 

3.60E-06 

3.46E-05 

1 

2.08E-04 

8.59E-05 

R21 

1.42E-04 

1.35E-05 

7.64E-05 

4.43E-05 

2.56E-06 

1.28E-05 

3.46E-05 

1 

3.67E-05 

R22 

1.11E-05 

1.44E-06 

7.48E-06 

1.11E-05 

2.01E-05 

2.56E-06 

3.61E-06 

9.16E-06 

1 

Table  6.  Non-diagonal  elements  in  the  T-matrix  normalized  to  diagonal  element  in  the  same  column  (Thin  shell  approximation) 
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c.  Evaluation  of  T -matrix  within  fluid  (at  radius  =  lm  and  2m) 

An  important  issue  is  to  ascertain  the  best  position  to  evaluate  the  T- 
matrix,  that  is,  on  the  shell  surface  or  at  some  distance  into  the  fluid.  To  investigate  this, 
the  T-matrix  diagonal  elements  were  computed  at  surfaces  within  the  fluid,  at  radii  of  1 
meter  and  2  meters,  respectively.  The  T-matrix  elements  detennined  by  COMSOL  are 
presented  in  Table  7,  with  the  errors  shown  in  the  last  two  columns: 


At  radius=lm 

I^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

0 

0 

P4 

-1.3698E-02 

1.1624E-01 

1.1704E-01 

9.6721E+01 

0.0078 

0.0002 

Rio 

-6.0761E-03 

7.773  IE-02 

7.7968E-02 

9.4470E+01 

-0.0916 

-0.0063 

R20 

-6.2362E-04 

-2.5007E-02 

2.5014E-02 

-9.1429E+01 

-0.2994 

0.0092 

At  radius=2m 

Rnm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

0 

0 

04 

-1.3699E-02 

1.1624E-01 

1.1704E-01 

9.6722E+01 

0.0077 

0.0009 

Rio 

-6.0763E-03 

7.773  IE-02 

7.7969E-02 

9.4470E+01 

-0.0908 

-0.0062 

R20 

-6.2379E-04 

-2.5012E-02 

2.5020E-02 

-9.1429E+01 

-0.2784 

0.0091 

Table  7.  COMSOL:  Diagonal  T-matrix  elements  at  radius  =  lm  and  2m  respectively  (fluid 
domain:  extra  fine  meshes  and  shell  surface:  extremely  fine  meshes) 


Comparing  to  the  T-matrix  evaluated  on  the  shell  surface  in  Table  4,  the 
magnitude  errors  for  R0o,  Rio  and  R20  were  0.001%,  -0.0449%  and  -0.1133% 
respectively.  This  has  obviously  shown  that  the  T-matrix  evaluated  on  the  shell  surface 
yields  a  more  accurate  result. 

C.  T-MATRIX  OF  A  THIN  SPHERICAL  SHELL  USING  3D  ELASTIC 
THEORY 

From  the  results  of  the  rigid  sphere  scattering  and  using  thin  shell  approximation 
for  the  shell,  we  have  developed  good  confidence  in  COMSOL  modeling.  The  next 
important  step  is  to  verify  whether  the  thin  shell  approximation  is  able  to  exhibit  the 
realistic  behavior  of  a  true  3D  elastic  shell.  This  involves  the  use  of  the  3D  elastic  theory 
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to  compute  the  T-matrix  of  a  fluid-loaded  thin  shell.  The  COMSOL  Multiphysics 
Acoustic-Structure  Interaction  (ACSI)  physics  mode  is  capable  of  performing  such 
simulation  to  achieve  accurate  values  of  the  T-matrix. 

1.  Analytical  T-matrix 

The  analytical  T-matrix  was  determined  based  on  the  3D  elastic  theory  and  the 
formulation  procedures  suggested  by  Hickling  [10].  The  T-matrix  diagonal  elements  are 
shown  in  Table  8.  The  last  two  columns  show  the  magnitude  error  percentage  in  the 
analytical  T-matrix  using  thin  shell  approximation,  relative  to  that  of  the  3D  elastic 
theory.  The  analytical  T-matrix  using  thin  shell  approximation  can  be  found  in  Table  3. 


^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Thin  shell 
Mag  Error 
(%) 

Thin  shell 
Phase  Error 
(deg) 

o 

o 

04 

-1.4191E-02 

1.1828E-01 

1.1913E-01 

9.6842E+01 

-1.76 

-0.12 

Rio 

-6.3614E-03 

7.9504E-02 

7.9758E-02 

9.4575E+01 

-2.15 

-0.10 

R20 

-6.21 12E-04 

-2.4915E-02 

2.4923E-02 

-9.1428E+01 

0.67 

-0.01 

Table  8.  Analytical  T-matrix  for  a  thin  shell  using  3D  elastic  theory 


From  the  magnitude  error  shown  in  Table  8,  it  is  obvious  that  the  thin  shell 
approximation  has  a  discrepancy  in  accuracy  of  up  to  2%,  relative  to  the  3D  elastic 
theory.  Hence,  it  is  important  to  note  that  computations  based  on  the  thin  shell 
approximation  theory  are  not  sufficiently  accurate  to  exhibit  the  real  behavior  of  a  3D 
element. 

2.  COMSOL  Results  for  3D  Elastic  Shell  Modeling 
a.  Diagonal  elements  of  T-matrix 

A  full  3D  elastic  steel  shell  model  was  created  in  the  COMSOL 
Multiphysics  Acoustic-Structure  Interaction  (ACSI)  physics  mode.  A  finer  mesh  was 
built  on  the  fluid  domain  and  a  coarser  mesh  for  the  solid  domain  -  the  entire  thickness  of 
the  shell.  The  element  size  is  as  follows: 
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Fluid  (finer  mesh): 


Maximum  element  size  =  0.275  m  (0.088A.) 

Shell  (coarser  mesh): 

Maximum  element  size  =  0.95  m  (95  x  thickness) 

The  mesh  of  the  model  consisted  of  164537  elements  and  is  computed 
with  318398  degrees  of  freedom  (DOF).  It  was  computed  over  an  average  time  frame  of 
7  minutes. 


The  diagonal  T-matrix  elements  employing  the  3D  elastic  theory, 
evaluated  on  the  shell  surface,  are  presented  in  Table  9.  The  last  two  columns  show  the 
magnitude  and  phase  errors  relative  to  the  analytical  T-matrix. 


^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

Roo 

-1.4187E-02 

1.1828E-01 

1.1912E-01 

9.6840E+01 

-0.0031 

-0.0019 

Ri,-i 

-6.3584E-03 

7.9489E-02 

7.9743E-02 

9.4573E+01 

-0.0185 

-0.0013 

Rio 

-6.3584E-03 

7.9489E-02 

7.9743E-02 

9.4573E+01 

-0.0188 

-0.0013 

Rn 

-6.3584E-03 

7.9489E-02 

7.9743E-02 

9.4573E+01 

-0.0185 

-0.0013 

R2,-2 

-6.1726E-04 

-2.4896E-02 

2.4904E-02 

-9.1420E+01 

-0.0766 

0.0078 

R2,-l 

-6.1720E-04 

-2.4895E-02 

2.4903E-02 

-9.1420E+01 

-0.0802 

0.0079 

R20 

-6.1723E-04 

-2.4896E-02 

2.4903E-02 

-9.1420E+01 

-0.0779 

0.0078 

R2I 

-6.1720E-04 

-2.4895E-02 

2.4903E-02 

-9.1420E+01 

-0.0802 

0.0079 

R22 

-6.1726E-04 

-2.4896E-02 

2.4904E-02 

-9.1420E+01 

-0.0766 

0.0078 

Table  9.  COMSOL:  Diagonal  T-matrix  elements  for  a  thin  shell  using  3D  elastic  theory 

(fluid  domain:  finer  meshes  and  shell:  coarser  meshes) 


An  overview  of  the  magnitude  errors  in  Table  9  shows  that  the  T-matrix 
detennined  by  COMSOL  using  the  3D  elastic  theory  is  in  excellent  agreement  with  the 
analytical  T-matrix.  It  has  also  shown  that  the  T-matrix  is  independent  of  the  azimuthal 
mode  number,  m. 

It  was  noted  that  the  mesh  of  the  3D  elastic  model  was  not  as  refined  as 


the  thin  shell  approximation  model.  Hence,  a  further  step  was  taken  to  examine  if  a  mesh 

28 


refinement  will  improve  the  results  of  the  3D  elastic  model.  A  test  case  on  R0o,  built  with 
finer  meshes  for  both  the  fluid  and  solid  domain  was  created.  The  computation  time  was 
2.5  hours  and  the  results  are  presented  in  Table  10. 


I^nm 

Real  part 

Imag  part 

Magnitude 

Phase  (deg) 

Mag 

Error 

(%) 

Phase 

Error 

(deg) 

o 

o 

04 

-1.4181E-02 

1.1825E-01 

1.19 10E-0 1 

9.6838E+01 

-0.0237 

-0.0032 

Table  10.  Improved  meshing  for  computing  Roo  using  3D  elastic  theory  (fluid:  finer  meshes 

and  solid:  finer  meshes) 

The  results  in  Table  10  show  that  the  magnitude  error  for  Roo  increases 
despite  the  mesh  refinement.  It  has  proven  that  mesh  refinement  may  not  necessarily  be 
the  improvement  factor  to  achieve  accuracy.  The  reason,  however,  is  unclear  at  this 
juncture  and  further  investigation  will  be  required. 

b.  Non-diagonal  elements  of  T-matrix 

The  non-diagonal  elements  relative  to  the  diagonal  elements  for  the  full 
3D  elastic  model  are  presented  in  Table  11.  The  largest  three  values  are  1.80E-03 
(monopole  to  quadrupole),  1.13E-03  (monopole  to  quadrupole)  and  1.11E-03 
(quadrupole  to  quadrupole).  These  values  will  not  cause  a  significant  leakage  in  the 
model.  Comparatively,  these  non-diagonal  values  are  lower  than  that  of  the  thin  shell 
approximation  model.  As  the  full  3D  elastic  shell  model  has  a  “real”  thickness  that 
enables  mesh  to  be  built  on,  the  computation  is  definitely  more  precise  than  the  model 
based  on  thin  shell  approximation. 
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Scattered  Pressure 


Incident  Pressure 


Roo 

Ri,i 

Rio 

Rn 

R2,-2 

R2,-l 

R20 

R21 

R22 

Roo 

1 

2.19E-05 

2.95E-05 

4.38E-05 

2.85E-06 

8.02E-06 

2.42E-06 

4.63E-05 

7.92E-05 

Ri-i 

1.99E-04 

1 

2.37E-05 

3.50E-05 

1.45E-05 

2.33E-05 

2.21E-05 

7.10E-05 

1.09E-04 

Rio 

1.32E-04 

1.17E-05 

1 

2.37E-05 

2.97E-06 

1.76E-05 

2.14E-05 

1.06E-04 

7.12E-05 

Rn 

9.95E-05 

8.95E-06 

1.17E-05 

1 

2.25E-06 

5.90E-06 

1.11E-05 

6.99E-05 

1.73E-04 

R2,-2 

1.80E-03 

7.69E-04 

3.90E-04 

6.10E-04 

1 

2.87E-04 

5.86E-04 

6.62E-04 

1.1  IE-03 

R2,-l 

1.13E-03 

3.68E-04 

4.87E-04 

3.81E-04 

7.17E-05 

1 

2.09E-04 

4.45E-04 

6.62E-04 

R20 

3.05E-05 

4.46E-05 

1.13E-04 

8.92E-05 

2.45E-05 

3.49E-05 

1 

2.09E-04 

5.86E-04 

R21 

1.90E-04 

3.16E-05 

8.10E-05 

1.23E-04 

4.59E-06 

1.24E-05 

3.49E-05 

1 

2.87E-04 

R22 

7.51E-05 

1.27E-05 

1.61E-05 

6.39E-05 

1.93E-06 

4.59E-06 

2.45E-05 

7.17E-05 

1 

Table  11.  Non-diagonal  elements  in  the  T-matrix  nonnalized  to  diagonal  element  in  the  same  column  (3D  Elastic  Theory) 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

The  T-matrix  of  a  fluid-loaded  thin  spherical  steel  shell  of  radius  a,  has  been 
calculated  using  the  COMSOL  Multiphysics  Finite-Element  Code  for  ka=  1  (where  k  is 
the  acoustic  wave  number).  In  order  to  observe  the  behavior  at  ka=  1,  the  thin  shell  radius 
was  chosen  to  be  0.5m,  with  a  thickness  of  1cm,  and  driven  at  a  frequency  of  474Hz.  The 
goal  of  attaining  a  result  accuracy  of  0.1%  relative  to  the  analytical  T-matrix  was 
achieved  in  a  model  employing  3D  elastic  finite  elements. 

The  modeling  process  was  structured  into  three  progressions.  The  starting  point 
was  to  compute  the  T-matrix  of  a  fluid-loaded  rigid  sphere  using  the  Acoustic-Shell 
Interaction  (ACSH)  physics  mode  to  examine  the  functionality  of  the  COMSOL 
Multiphysics  code.  Secondly,  using  the  same  physics  mode  ACSH,  the  T-matrix  of  a 
fluid-loaded  thin  spherical  shell  was  determined  by  employing  the  thin  shell 
approximation  theory.  Lastly,  the  T-matrix  of  a  true  elastic  shell  was  determined  by 
adopting  the  3D  elastic  theory  using  the  Acoustic-Structure  Interaction  (ACSI)  physics 
mode. 

The  T-matrix  elements  for  the  rigid  sphere,  evaluated  on  the  sphere  surface,  were 
found  to  be  marginally  satisfactory  relative  to  the  analytical  values.  The  quadrupole,  R20 
was  observed  to  be  the  worst  case  with  the  highest  magnitude  error.  The  magnitude  error 
with  respect  to  the  analytical  values  showed  an  improvement  from  -1.032%  to  -0.335% 
when  a  mesh  refinement  of  the  fluid  and  solid  domain  was  introduced. 

In  the  evaluation  of  the  T-matrix  on  the  shell  surface  using  thin  shell 
approximation  theory,  the  results  appeared  to  be  in  better  agreement  than  those  of  the 
rigid  sphere.  The  T-matrix  magnitude  error  were  ranged  from  -0.28%  to  0.001%,  relative 
to  the  analytical  values  for  a  thin  shell.  The  non-diagonal  elements  had  shown  that  there 
was  no  significant  “leakage”  in  the  model.  Also,  the  T-matrix  was  evaluated  at  two  radii 
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(1  meter  and  2meters)  within  the  fluid.  The  results  determined  were  not  as  favorable  as 
those  evaluated  on  the  shell  surface. 

In  modeling  the  thin  shell  as  a  true  3D  elastic  element,  the  simulation  based  on  the 
3D  elastic  theory  was  able  to  achieve  the  goal  of  attaining  an  accuracy  of  at  least  0.1%. 
The  model  was  computed  using  the  Acoustic-Structure  Interaction  (ACSI)  physics  mode 
and  the  T-matrix  was  evaluated  on  the  shell  surface.  The  resulting  average  magnitude 
error  was  less  than  0.08%,  relative  to  the  analytical  values.  Also,  the  non-diagonal  values 
were  not  significant. 

The  COMSOL  Multiphysics  Finite-Element  Code  has  proven  its  ability  to  model 
the  scattering  of  acoustic  waves  from  a  thin  shell.  Comparing  the  T-matrix  results 
obtained  by  ATILA  in  the  previous  attempt,  COMSOL  is  superior  in  tenns  of  efficiency 
and  accuracy. 

B.  FUTURE  WORK 

In  spite  of  the  satisfying  results  obtained  in  COMSOL  Multiphysics  Finite- 
Element  Code,  there  are  three  areas  that  can  be  explored  to  further  improve  the  results 
accuracy. 

1.  External  determination  of  scattering  coefficients 

It  is  unclear  about  the  preciseness  of  the  post  evaluation  performed  by  COMSOL 
which  can  possibly  introduce  small  numerical  errors.  It  can,  therefore,  be  investigated  by 
performing  an  external  determination  of  the  scattering  coefficients  by  an  alternate 
program  such  as  MATLAB.  In  this  case,  the  nodal  values  obtained  in  COMSOL  will 
have  to  be  exported  to  MATLAB  in  which  a  singular  value  decomposition  using 
Equation  (28)  can  be  perfonned  by  the  least  square  curve-fitting  procedures. 

2.  Perform  swept  mesh  in  radial  direction 

Since  a  refined  mesh  may  not  be  necessary  for  the  entire  domain,  for  instance,  the 

fluid  domain,  the  swept  mesh  technique  can  be  explored.  The  swept  mesh  allows 

different  refinement  on  a  surface.  With  this  option,  the  user  will  be  able  to  specify  the 
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essential  size  of  the  mesh  and  sweep  radially  within  the  domain.  In  this  case,  the 
computing  resource  can  be  preserved. 

3.  Perfectly  Matched  Layer  (PML) 

The  absorption  of  outgoing  waves  can  be  investigated  by  replacing  the  radiation 
boundary  condition  at  the  fluid  outer  radius  by  a  perfectly  matched  layer  (PML)  to 
observe  the  difference  in  wave  scattering  amplitudes.  A  PML  may  perfonn  better  in 
tenns  of  absorbing  outgoing  radiation  when  the  scatterer  is  not  spherical.  This  remains  to 
be  investigated. 
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